clear all;
clc;
folder='sigma_29Apr08';

%
tmax=1000;
T=1:tmax;

sig_0_025 = dlmread(strcat(folder,'\sig_0_025.dat'));
sig_0_025_up = dlmread(strcat(folder,'\sig_up_0_025.dat'));
sig_0_025_down = dlmread(strcat(folder,'\sig_down_0_025.dat'));


sig_0_05 = dlmread(strcat(folder,'\sig_0_05.dat'));
sig_0_05_up = dlmread(strcat(folder,'\sig_up_0_05.dat'));
sig_0_05_down = dlmread(strcat(folder,'\sig_down_0_05.dat'));


sig_0_1 = dlmread(strcat(folder,'\sig_0_1.dat'));
sig_0_1_up = dlmread(strcat(folder,'\sig_up_0_1.dat'));
sig_0_1_down = dlmread(strcat(folder,'\sig_down_0_1.dat'));


sig_0_2 = dlmread(strcat(folder,'\sig_0_2.dat'));
sig_0_2_up = dlmread(strcat(folder,'\sig_up_0_2.dat'));
sig_0_2_down = dlmread(strcat(folder,'\sig_down_0_2.dat'));


figure(1)
plot(T,sig_0_025,T,sig_0_05,T,sig_0_1,T,sig_0_2);
xlabel('T:1000');
ylabel('\sigma_\epsilon(t)');
title('\sigma_\epsilon(t) vs T for \epsilon = 0.025, 0.05 (200 iterations), 0.1 and 0.5 for 1000 iterations');
grid on;

% figure(2)
% plot(T,sig_0_025,T,sig_0_025_up,T,sig_0_025_down);
% xlabel('T:1000');
% ylabel('\sigma_\epsilon(t)');
% title('\sigma_\epsilon(t) vs T for \epsilon = 0.025 for charge and up and down spins');
% grid on;

% figure(2)
% plot(T,sig_0_025);
% 
% figure(3)
% plot(T,sig_0_05);
% 
% figure(4)
% plot(T,sig_0_1);
% 
% figure(5)
% plot(T,sig_0_2);
